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Temperature  Dependence  of  Low  Energy  Carrier 
Dynamics  of  Silicon  by  Terahertz  Time  Domain 

Spectroscopy 

Shigeki  Nashima,  Osamu  Morikawa,  Kazue  Takata,  Masanori  Hangyo, 


Abstract  - The  complex  conductivity  of  moderately  P-doped 
silicon  wafers  (1.1±0.2  Clem)  at  various  temperatures  has 
been  deduced  by  using  a terahertz  time-domain  (THz-TDS) 
spectroscopy.  The  characteristic  dispersion  of  the  complex 
conductivity  is  observed  in  the  THz  region.  Temperature 
dependence  of  the  complex  conductivity  changes 
dramatically  around  —120  K,  which  is  interpreted  in  terms 
of  enhancement  of  mobility  and  the  freezing  of  the  carrier. 
The  experimental  data  at  low  temperatures  slightly  deviates 
from  the  simple  Drude  model,  which  causes  the 
underestimation  of  the  mobility  and  overestimation  of  the 
carrier  density. 

I.  Introduction 

It  is  well  known  that  the  optical  properties  of 
semiconductors  from  microwave  to  infrared  regions  are 
greatly  influenced  by  the  dynamics  of  free  carriers  of 
semiconductors!  l]-[6].  Dispersion  of  the  complex 
conductivity  (crr/cr2)  depends  on  the  carrier  density  and 
scattering  time,  which  is  explained  by  the  simple  Drude 
model.  Such  a model  conversely  provides  a simple 
method  of  calculating  optical  constants  in  the  far  infrared 
region  from  the  carrier  density  and  mobility.  However, 
the  optical  properties  in  this  region  at  low  temperatures 
have  not  been  discussed  so  much  before.  The  terahertz 
time-domain  spectroscopy  (THz-TDS)  is  used  to  measure 
the  complex  conductivity  of  samples  at  THz  frequencies 
without  the  need  of  the  Kramers-Kronig 
transformation^].  In  this  letter,  the  complex  conductivity 
of  a moderately  doped  silicon  wafer  is  deduced  at  various 
temperatures  by  using  the  THz-TDS.  We  obtain  the 
carrier  density  and  mobility  of  the  sample  at  low 
temperatures  and  compare  them  with  the  results  obtained 
from  the  Hall  measurement  and  four-point  contact 
resistivity. 

II.  Experimental 


An  tt-type  silicon  wafer  doped  with  phosphorus  was  used 
as  a sample  in  this  study.  The  resistivity  at  room 
temperature  and  the  thickness  of  the  sample  were 
1.1  (±0.2)  Qcm  and  400  pm,  respectively.  A conventional 
THz-TDS  system  was  used  wherein  the  pump  and  the 
trigger  beams  were  derived  from  a mode-locked 
Tirsapphire  laser  which  produces  optical  pulses  with  a 
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width  of  50  fs  and  a wavelength  of  790  nm  at  a repetition 
a rate  of  82  MHz[4][5].  The  sample  was  mounted  on  a 
sample  holder  of  a closed-cycle  He  refrigerator  specially 
designed  to  avoid  mechanical  vibrations  and  was  inserted 
between  a set  of  paraboloidal  mirrors.  In  this 
measurement,  a 5-pm  gap  bowtie-type  photoconductive 
antenna  was  used  as  a THz  emitter  and  detector.  The 
frequency  range  of  this  system  was  limited  by  the 
absorption  property  of  a 3-mm-thick  fused  quartz 
windows  of  the  He  refrigerator  and  was  estimated  to  be 
less  than  1.5  THz. 


III.  Results  and  Discussion 

We  calculate  the  complex  refractive  index  n from 
experimental  transmittance  and  phase  shift  obtained  by 
the  Fourier  transformation  of  the  THz  pulse  waveforms 
before  and  after  inserting  the  sample.  The  complex 
conductivity  is  obtained  from  the  complex  refractive 
index  from  the  relation, 

n2  -e(co)-£^  +icr(a>)/ a>£0,  (1) 

where  £(£]-ie2)  is  the  complex  dielectric  constant,  % is 
the  dielectric  constant  of  undoped  silicon  and  £0  is  the 
free-space  permittivity.  We  used  the  expression 

£sj=(3.415)2[l].  Frequency  dependence  of  the  complex 
conductivity  of  the  sample  at  300  K,  120  K,  70  K,  and  16 
K is  shown  in  Figs.  2(a)-(d),  respectively.  At  room 
temperature  the  value  of  cr]  decreases  monotonically  with 
increasing  frequency.  The  value  of  a2  at  room 
temperature  increases  from  zero  with  increasing 
frequency  until  ~0.6  THz  and  decreases  monotonically 
with  increasing  frequency  further.  At  0.6  THz,  a 
crossover  of  the  two  curves  is  observed.  Frequency 
dependence  of  the  complex  conductivity  becomes 
pronounced  with  decreasing  temperature  until  -120  K. 
The  value  of  cr,  in  the  low  frequency  region  increases  and 
that  in  the  high  frequency  region  decreases  with 
decreasing  temperature.  The  temperature  dependence  of 
cr2  also  becomes  pronounced  with  decreasing  temperature. 
The  crossover  point  of  the  two  curves  shifts  to  the  lower 
frequency  region  with  decreasing  temperature.  The 
dispersion  of  the  complex  conductivity  decreases  rapidly 
with  decreasing  temperature  below  120  K and  the  values 
of  <7j  and  <x2  at  16  K are  almost  zero  over  the  entire 
frequency  region. 
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Solid  lines  in  Fig.  1 show  the  fitting  curves  to  cr,  and  <j2 
at  each  temperature  using  the  simple  Drude  model. 
According  to  the  simple  Drude  model,  in  which  the 
collision  damping  is  independent  of  the  carrier  energy, 
the  frequency-dependent  complex  conductivity  in  SI  units 
is  given  by 


ex  = i£0oS~  !{co  + if) , (2) 


where  r-Mr  is  the  carrier  damping  rate  and  r is  the 
average  collision  time.  The  plasma  frequency  &>p  is 
defined  by  co2~Nce2/£0m\  where  Nc  is  the  number  density 
of  carriers,  e is  the  electronic  charge  and  m is  the 
effective  carrier  mass.  The  expression  m*=0.26  w0  is  also 
used,  where  m0  is  the  free-electron  mass[7].  Our  fitting 
procedure  effectively  optimizes  the  values  of  Nc  and  r, 
from  which  the  mobility  p given  by  p=er/w*  is 
determined.  The  two  parameters  are  determined  by  fitting 
the  complex  conductivity  data  between  0.2  and  1 THz. 
Figure  2 (a)  shows  the  temperature  dependence  of  the 
mobility  obtained  from  the  fitting  (open  circles).  Solid 
and  dashed  lines  show  the  data  obtained  from  previous 
Hall  measurements  on  two  samples  with  smaller  and 
larger  resistivities  (p  = 0.66  Qcm,  4.76  Qcm)  than  the 
present  sample,  respectively[8j.  Error  bars  show  the 
uncertainties  associated  with  temperature.  At  room 
temperature  the  value  of  p and  p (=l/a(0))  is  estimated  to 
be  1560  cm2/Vs  and  1.06  Qcm  for  the  present  sample, 
respectively,  which  agree  on  the  whole  with  that  obtained 
by  Hall  measurement  and  four-point  contact  resistivity 
measurement^].  The  value  of p increases  dramatically 
with  decreasing  temperature  until  150  K.  The  increasing 
rate  gradually  decreases  with  decreasing  temperature. 
However,  the  value  of  p obtained  by  the  simple  Drude  fit 
at  various  temperatures  is  a little  smaller  than  that  of  the 
silicon  with  lower  resistivity  obtained  from  the  Hall 
measurement. 

Figure  2 (b)  displays  the  temperature  dependence  of  the 
carrier  density  (open  circles)  of  the  sample  obtained  from 
the  simple  Drude  fit.  According  to  the  Fermi-Dirac 
distribution,  the  carrier  density  of  moderately  doped 
silicon  at  various  temperatures  satisfies  the  following 
equation, 


Nn-N. 
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where  N0  is  the  donor  density,  Nu  is  the  number  of  the 
valley  at  the  X point,  mn  is  the  effective  mass,  kB  is  the 
Boltzman  constant,  T is  the  temperature,  h is  the  Planck 
constant,  and  AED  is  the  donor  ionization  energy  of  the 
phosphorus  dopant.  The  value  of  Nu  is  6.  The  values, 
Ad=3.4x1015  cm'3  (solid  lines)  and  5xl015  cm'3  (dashed 
lines),  are  used  corresponding  to  1.1 +0.2  Qcm  and  1.1- 
0.2  Qcm  (the  positional  fluctuation  of  the  resistivity  of 
the  present  sample)[10].  The  expression 
m=(m2m])m:=032m0  is  also  used,  where  mx  and  mx  are 
longitudinal  and  transverse  masses  at  the  X point, 
respectively.  At  room  temperature,  the  experimentally 


Fig.  1.  Temperature  dependence  of  the  complex 
conductivity  a\  (open  circles)  and  a2  (solid  circles)  at 
(a)  300  K,  (b)  120  K,  (c)  70  K,  and  (d)  16  K.  Solid 
lines  show  the  simole  Drude  model. 


Temperature  (K) 

(a) 


(b) 

Fig.  2.  Temperature  dependence  of  (a)  mobility  and  (b) 
carrier  density  deduced  from  the  simple  Drude  fit  to  the 
complex  conductivity  obtained  from  the  THz-TDS 
measurement.  The  solid  and  dashed  lines  in  (a)  show 
the  results  from  Hall  measurements  for  silicon  with 
resistivity  of  0.66  Qcm  and  4.76  Qcm,  respectively.  The 
solid  and  dashed  lines  in  (b)  show  the  carrier  density 
calculated  using  Eq.  (3)  assuming  AD=3.4xl015  cm*3 
and  5x1 015  cm'3  as  a function  of  temperature, 
respectively. 

obtained  Nc  is  found  to  be  3. 7xl0r 5 cm'3  for  the  present 
sample.  The  experimentally  obtained  value  of  Nc 
increases  slightly  with  decreasing  temperature  until  150 
K and  rapidly  drops  to  zero.  Below  30  K,  the  accurate 
estimation  of  the  carrier  density  from  the  experimental 
results  is  difficult  because  of  small  frequency  dependence 
of  complex  conductivity  compared  to  the  noise  level. 

At  room  temperature,  the  results  obtained  from  THz-TDS 
agree  well  with  the  simple  Drude  model.  However,  the 
deviation  from  the  simple  Drude  model  is  observed  at 


low  temperatures  as  shown  in  Figs.  l(a)-(d).  The  value 
of  ctj  at  low  frequency  region  obtained  from  THz-TDS  is 
smaller  than  that  obtained  from  the  simple  Drude  model, 
which  is  similar  to  the  results  reported  by  Exter  et  al[ 4]. 
Deviation  from  the  simple  Drude  model  becomes  more 
pronounced  with  decreasing  temperature.  We  calculate 
ju(a))  and  Nz(co)  from  the  experimental  complex  dielectric 
constant  at  each  frequency  at  various  temperatures  by  the 
following  equations[2]. 

/<(*)= (4) 

m (0£2 

Arc(«)=^ho^(1+«2^)}’  (5) 

e z 

The  dispersion  of  jj{co)  and  Nc(co)  at  various  temperatures 
are  plotted  in  Figs.  3 (a)  and  (b),  respectively.  At  room 
temperature,  the  value  of  jj(co)  is  constant  over  the  entire 
frequency  region,  which  indicates  that  the  experimental 
results  agree  well  with  the  simple  Drude  model.  However, 
the  values  of  jj{co)  at  low  temperatures  decrease  with 
increasing  frequency,  which  means  disagreement  with  the 
simple  Drude  model.  This  dispersion  is  also  observed  at 
90  K in  the  millimeter  wave  region  and  is  explained  by 
the  calculation  based  on  the  Boltzman  transport 
equation[3].  The  dispersion  of  /j{(o ) becomes  more 
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Fig.  3.  The  dispersion  of  (a)  fj.  and  (b)  Nc  obtained  from 
the  experimental  dispersion  of  and  £2  at  300  K (open 
circles),  120  K (solid  circles),  70K  (open  squares),  and 
16  K (solid  squares). 


pronounced  with  decreasing  temperature.  Below  50  K 
jj(co)  is  difficult  to  estimate  because  of  small  frequency 
dependence  compared  to  the  noise  level.  These  results 
indicate  that  jj(cd=  0)  is  larger  than  jj  obtained  from  the 
simple  Drude  fitting  in  the  entire  frequency  region. 
Although  the  value  of  Nc(co)  at  room  temperature  is  also 
constant  as  well  as  //(&>),  the  change  of  Nc(&)  with 
frequency  is  observed  at  low  temperatures.  The  value  of 
Nc(co)  slightly  increases  with  increasing  frequency,  which 
predicts  that  Nc  at  low  temperatures  obtained  from  the 
simple  Drude  model  estimates  a slightly  larger  than  that 
obtained  from  Nc(o)  = 0). 


IV.  Conclusion 

The  THz-TDS  system  is  used  to  measure  the  complex 
conductivity  of  n- type  silicon  wafer  at  temperatures  from 
16  K to  room  temperature.  The  measured  complex 
conductivity  is  well  fitted  by  the  simple  Drude  model  at 
room  temperature  and  the  experimentally  obtained  carrier 
density  and  mobility  agree  on  the  whole  with  Hall 
measurement  and  four-point  contact  resistivity.  This 
result  indicates  that  the  THz-TDS  measurement  is  an 
effective  method  for  electric  noncontact  characterization 
of  doped  silicon.  The  change  of  the  temperature 
dependence  of  the  complex  conductivity  around  -120  K 
is  observed  and  slight  deviation  from  the  simple  Drude 
model  is  found  at  low  temperatures.  This  deviation 
causes  the  calculation  of  the  mobility  and  carrier  density 
using  the  simple  Drude  model  in  the  entire  frequency 
region  to  be  slightly  underestimated  and  overestimated, 
respectively. 
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